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Validity of Stokes-Einstein Relation in Soft Colloids up to the Glass Transition 
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We investigate the dynamics of kinetically frozen block copolymer micelles of different softness 
across a wide range of particle concentrations, from the fluid to the onset of glassy behavior, through 
a combination of rheology, dynamic light scattering and pulsed held gradient NMR spectroscopy. We 
additionally perform Brownian dynamics simulations based on an ultrasoft coarse-grained potential, 
which are found to be in quantitative agreement with experiments, capturing even the very details of 
dynamic structure factors S{Q,t) on approaching the glass transition. We provide evidence that for 
these systems the Stokes-Einstein relation holds up to the glass transition; given that it is violated 
for dense suspensions of hard colloids, our Hndings suggest that its validity is an intriguing signature 
of ultrasoft interactions. 

PACS numbers: 82.70.-y, 61.20.Gy, 64.70.km, 82.70.Dd 


The microscopic origin of vitrification as a system un¬ 
dergoes quenching either by a sudden decrease irr tem¬ 
perature T, as in supercooled liquids [1-3], or by a fast 
increase of the volume fraction /, as in colloidal glasses 
[4-6], is still a subject of intense discussion. Colloidal 
systems have been playing an important role in unveiling 
the microscopic aspects of the transition due to the ca¬ 
pability of tuning specific interactioirs and in developing 
experimental and simulation techniques that permit in¬ 
vestigations both at the single-particle and the collective 
levels. Accordingly, it is highly desirable to have at one’s 
disposal a versatile and well-controlled experimental sys¬ 
tem, which allows for a systematic tuning of softness at 
the individual particle level. In addition to fundamental 
interest, tunable rheological behavior of soft colloids [7] is 
very important for tailoring material properties with rel¬ 
evance to technological applications in oil and medical in¬ 
dustries [8] . Recently, kinetically-frozen block-copolymer 
micelles have emerged as easy-to-synthesize and precisely 
tunable model systems for soft-colloid suspensions [9-14]. 
Since the softness and morphology of the micelles can be 
changed by several factors, which include the aggregation 
number, the solvent composition [9], the solvophobic-to- 
solvophilic block ratio or the block length, these micelles 
allow to systematically bridge the physics of linear poly¬ 
mers to that of colloidal hard spheres. 

One of the fingerprints of a (metastable) liquid ap¬ 
proaching the glass transition is the breakdown of the 
Stokes-Einstein (SE) relation [15]: for glass-forming sys¬ 
tems the product of the macroscopic (zero-shear) viscos¬ 
ity r] and the mesoscopic (long-time) self-diffusion coef¬ 
ficient Dg is not constant anymore as the system slows 
down. The origin of this behavior has been attributed to 
the emergence of dynamic heterogeneities, provided by 


the presence of distinguishable populations of fast and 
slow particles [16-19]. The breakdown of the SE relation 
close to the glass transition is a generic feature shared by 
atomic [20] and molecular [21-24] glass formers, polymers 
[25] and metallic glasses [26]. Interestingly, only a lim¬ 
ited number of experimental studies on colloidal systems 
showing deviations from SE behavior is available, in¬ 
cluding hard spheres [27] and Laponite suspensions [28]. 
Much more abundant in the literature are numerical re¬ 
sults, which have allowed to grasp microscopic insights 
on the nature of the SE breakdown and have shown that 
both the onset and the strength of the deviations depend 
on the specific interaction potential. Thus, while atomic 
and molecular glass-formers inevitably display large de¬ 
viations upon lowering temperature, colloidal systems, 
being amenable of a systematic and controlled variation 
of effective interactions, offer the ideal playground for as¬ 
sessing whether a more extended range of validity of SE 
exists for specific systems. 

In hard spheres suspensions, systematic SE violations 
were found both in experiments [27] and in simulations 
[29]; moreover, for short-range attractive colloids an en¬ 
hanced breakdown was observed, a fact that has been 
rationalized in terms of pronounced dynamical hetero¬ 
geneities [29-31]. In contrast, for very soft potentials like 
the Gaussian core model, no significant violation of SE 
was observed [32, 33], suggesting that such models are 
“mean-field”, in the sense that they suppress dynamic 
heterogeneities due to the presence of long-range interac¬ 
tions. These observations are also very relevant to the re¬ 
cent discussion of whether there exists special limits, such 
as high dimensions, where a mean-field description of the 
glass transition may hold [32, 34-36]. Indeed numerical 
simulations of hard hyperspheres have also shown a sys- 
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TABLE I: Molecular weight [kg/mol], block ratio m\n of 
PBO to PEO and PEP to PEO repeat units, polydispersity 

o 

p, aggregation number A^agg, and sizes Ru [A] (hydrodynamic 

o 

radius) and (Tint [A] (micelle interaction diameter). 


Diblock 

copolymer 

Mu, 

PBO 

Mu, 

PEP 

Mu, 

PEO 

Block 

ratio 


-Yagg 

Rh 

(Tint 

hPBOlO-dPEOSO 

10.5 

- 

52.5 

1:8 

1.02 

500 

727 

899 

hPEPl-hPEO20 

- 

1.1 

21.9 

1:33 

1.04 

120 

334 

380 


“Afuj/Mn polydispersity determined by ^H-NMR and SEC. 


tematic reduction of SE violation with increasing dimen¬ 
sionality [37, 38]. The validity of SE relation up to the 
onset of glassy arrest may have important consequences 
also for the experimental community, as it would provide 
a test-case for the use of advanced microrheological tech¬ 
niques also close to arrest [39]. Thus, the identification 
of experimental systems which strictly satisfy SE even in 
the glassy regime have a broad interest in the scientific 
community. 

In this Letter, we report a comprehensive study of 
frozen block copolymer micelles featuring tunable soft¬ 
ness [9, 11, 12]. Experimental results from rheology mea¬ 
surements, dynamic light scattering (DLS), and pulsed 
field gradient NMR spectroscopy (PFG-NMR) are com¬ 
pared to those obtained from Brownian dynamics sim¬ 
ulations (BD) based on a coarse-grained description of 
the system. We find excellent agreement between experi¬ 
ments and theoretical description for the dynamic struc¬ 
ture factor S{Q, t) over a broad range of volume fractions 
approaching the glass transition value 4>g, as well as for 
the density dependence of the transport coefficients, in¬ 
cluding the self-diffusion coefficient Dg , the shear viscos¬ 
ity rj and the collective relaxation time r. Most strik¬ 
ingly, we find that the SE relation remains valid in the 
whole investigated range, up to the glass transition, for 
two different values of the softness parameter, suggesting 
that star-like micelles may be viewed as an experimental 
realization of mean-field glasses, where dynamic hetero¬ 
geneities are strongly suppressed. 

We investigate aqueous solutions of star-like micelles 
with two representative aggregation numbers cover¬ 
ing the regime from typical ultra-soft (A^agg = 120) 
to moderately soft (iVagg = 500) colloidal particles. 
For the first system, the amphiphilic block copolymer 
poly(ethylene-alt-propylene)-poly(ethylene oxide), PEP¬ 
PED, was employed, while the second was formed by 
poly(butyleneoxide)-poly(ethyleneoxide), PBO-PEO [11, 
14]. We covered a broad range of polymer volume frac¬ 
tion (0.06% < (f) < 4%) from very dilute suspensions to 
well above the glass transition for both systems. The 
samples’ characteristics are summarized in Table I. 

To describe the system theoretically, we employ a 
coarse-grained model, originally developed for regular 
star polymers in good solvent conditions [40, 41], which 
was recently shown to be able to describe the structure 


and phase diagram of star-like micelles [11]. In this 
framework, micelles are coarse-grained as soft particles 
with an effective, repulsive potential 14ff(r) acting be¬ 
tween their centers. The latter features a logarithmic 
dependence at short distances (r < (Tint) and a Yukawa 
decay at large ones (r > (Tint)- The so-called corona di¬ 
ameter (Tint, roughly equal to the hydrodynamic radius, 
sets the interaction range while the aggregation number 
A^agg, analogous to the star functionality / [40, 41], tunes 
the softness of the repulsion. This model allows to de¬ 
scribe the phase behavior of the micellar suspension on 
the space spanned by A^agg and the micellar concentration 
0TH = (7’'/6)pmtTinti where pm = NmlV is the density of 
a system containing Nm micelles in the volume V. While 
A(j,gg is independently determined by SANS form factor 
analysis in dilute solution [9], pm can be unambiguously 
determined in experiment in terms of weighted samples. 
The correspondence between pm and ^th is established 
via the experimental hydrodynamic radius Rh, which is 
used to determine (Tint, as previously described in equi¬ 
librium studies of these systems [11]. To study the long¬ 
time dynamics of the model, we simulate = 2000 
spherical particles of unit mass to, immersed in a cubic 
box of fixed volume Vbox, interacting through the coarse¬ 
grained potential 14ff(r). To avoid crystallization at high 
densities, particle sizes are drawn from a Gaussian distri¬ 
bution of standard deviation 10% and average diameter 
(Tint, which sets the unit of length. Simulations are per¬ 
formed at fixed temperature, setting fcsT = 1 with ks 
the Boltzmann constant, for A^agg = 120 and A^agg = 500 
at various volume fractions (pTH controlled by changing 
Vbox- To mimic the solvent effectively, we use Brown¬ 
ian dynamics (BD) simulations with the bare, short-time 
diffusion coefficient fixed to Dq = 0-001 [42]. The inte¬ 
gration time step is ch osen as At = 0.01, where time is 
measured in units of ymaf^^JJkBT). 



FIG. 1: Theoretical phase diagram [41, 43] and correspond¬ 
ing experimental path studied in this work. Open and closed 
symbols indicate investigated fluid and glassy samples, re¬ 
spectively. Triangles mark the glass transition obtained by 
extrapolating diffusivity data to zero through a MCT-power- 
law. The cartoons illustrate the different softness of the mi¬ 
celles at the two studied Nagg. 
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Figure 1 displays the phase diagram of the system in 
the (0TH, -/Vagg)-plane [41], comparing the experimen¬ 
tal state points investigated in this work with theoreti¬ 
cal results for the fluid-solid [41] and glass (MCT) [43] 
boundaries. The estimates of the polymer volume frac¬ 
tions at the glass transition, determined by rheology, are 
= (3.476 ±0.03)% and (2.635 ±0.05)% for IVagg = 120 
and 500 respectively, corresponding to micellar packing 
fractions ((iTH.g — 0.236±0.006 and 0.230±0.010, in good 
agreement with MCT predictions [43]. 

Let S{Q,t) = exp[-iQ • (r,(<) - rg(0))]) be 

the dynamical structure factor of a system of N parti¬ 
cles, rk{t) denoting the position of particle k at time t, 
and S{Q) the equal-times, static structure factor. A de¬ 
tailed comparison between the ratios S{Q,t)/S{Q) mea¬ 
sured by DLS and those calculated from BD simulations 
is reported in Fig. 2 for different volume fractions at 
fixed wavevectors Qa-mt — 1-68 and 1.07 for A^agg = 500 
and 120, respectively. To superimpose simulation data 
onto the experimental timescale, the microscopic dynam¬ 
ics is adjusted through an arbitrary shift of the time- 
axis, which depends on A^agg but is the same for all vol¬ 
ume fractions. As ^th increases, the typical pattern of 
glass-forming systems emerges in S{Q,t)\ a two-step re¬ 
laxation develops with a growing plateau at intermedi¬ 
ate times followed by a long-time final relaxation, which 
grows on approaching the glass transition. The experi¬ 
mentally measured correlators are found to be somehow 
intermediate (e.g., in the value of the plateau) between 



FIG. 2: Time dependence of normalized dynamic structure 
factors from DLS (symbols) compared with self correlators 
calculated from BD simulation (solid lines). Top panel: Re¬ 
sults for QtJint cs 1.68 and Akgg = 500. As a reference, also 
the corresponding collective correlator (dashed line) calcu¬ 
lated from BD at the same Quint is reported for the largest 
(j>TH- Bottom panel: Same as before for Quint — 1.07 and 
A^agg = 120. Simulation data were multiplied by a constant 
factor of 11 fis and 2 fis for 77agg = 500 and 120, respectively, 
to match the experimental microscopic relaxation. 



FIG. 3: Wavevector dependence of the normalized dynamic 
structure factors from DLS (symbols) and collective density 
correlators from BD simulation (solid lines) for (j)TH = 0.17 
and A^agg = 120 (left panel), and for (j>TH = 0.2 and 77agg = 
500 (right panel). 


the self and collective numerical ones (see Supplementary 
Information [44]). However, the long-time relaxation is 
well described by both the self and collective numeri¬ 
cal data, capturing quantitavely the growth of the relax¬ 
ation time with increasing (()th a-nd the shape (i.e., the 
stretching exponent) of the correlators. We stress that 
apart from the adjustment of the microscopic time, there 
is no fit parameter in the comparison, yielding an al¬ 
most quantitative description within experimental error 
for both studied values of N^gg- 

Additional evidence for the accuracy of the coarse¬ 
grained interaction to describe the equilibrium dynamics 
is provided by the comparison of the S{Q,t)/S{Q) from 
experiments and simulations as a function of the scatter¬ 
ing vector Q, which is shown in Fig. 3 at constant volume 
fraction. Clearly, the Q-dependence of the slowest relax¬ 
ation process in terms of characteristic time and shape 
(stretching exponent), as well as the Q-dependence of the 
intermediate non-ergodicity plateau are well described by 
BD simulations for both values of A^agg, covering a length 
scale variation between 0.5 < Qaint ^ 2.2; as a reference, 
close to the glass transition the nearest-neighbor peak of 
S{Q) is found at Qa-mt — 5. 

So far we have shown that the coarse-grained poten¬ 
tial [40] can describe both the structure [11, 13] and the 
dynamics, in a semi-quantitative way, of frozen block 
copolymer micelles. Building on this, we now provide 
evidence of the validity of SE relation in these systems 
by combining direct experimental measurements and in¬ 
formation that can be extracted from the theoretical de¬ 
scription of the dynamical correlators. To this end, we 
have investigated the concentration dependence of the 
long-time, self-diffusion coefficient Ds{(j)) by means of 
both DLS and PFG-NMR, and also the zero-shear viscos¬ 
ity r]{cj)) through rheological measurements. In Fig. 4 the 
reduced viscosity rj[(j))/riQ and the inverse of the reduced 
self-diffusion coefficient DQ/Ds{(j>) are reported as a func¬ 
tion of approaching the glass transition for both ex¬ 
perimental systems with A^agg = 120 and 500. Here rjo 
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FIG. 4: Inverse of reduced self-diffusion coefficient (Dq/Ds) 
from DLS, PFG-NMR and BD simulations, reduced zero- 
shear viscosity (7?/?7o) from rheology, and reduced relaxation 
time (r*/ro) calculated at Quint = 5.0 from BD simulations, 
are shown together as a function of volume fraction (pTH for 
Na,gg = 500 (black) and A^agg = 120 (red). The experimental 
glass transition for both A^agg are indicated by the vertical 
solid lines. The corresponding MGT glass lines, obtained by 
power law fits of the numerical data, are indicated by vertical 
dashed lines. Inset: DsT] from experiments and DsT* from 
simulations as function of the distance to the glass transition 
e. Gollective relaxation times r* were calculated by fitting 
the correlators at the peak of S{Q) to stretched exponentials. 
Hard-sphere data for DsT) were obtained by MD simulations 
as explained in Ref. [29]. 

represents the solvent viscosity, whereas Dq is the diffu¬ 
sion coefficient of the micelles at infinite dilution. For 
comparison, in Fig. 4 we also report results for the long¬ 
time diffusion coefficient Ds and the collective structural 
relaxation time r extracted from BD simulations. Ex¬ 
ploiting the fact that the viscosity is proportional to the 
relaxation time, we present data for the latter at fixed 
wavevector Quint = 5 and normalized by the correspond¬ 
ing time at infinite dilution, tq. Excellent overlap (within 
experimental error) between measured data and results 
from BD simulation based on the coarse-grained model is 
found without the introduction of any adjustable param¬ 
eter. Furthermore, we establish that up to the largest 
volume fractions in Fig. 4, the self-diffusion coefficient, 
the viscosity and the relaxation time for each of the two 
studied systems follow each other very closely, approach¬ 
ing the glass transition in a similar way. 

These findings suggest that SE holds for such soft block 
copolymer micelles with star-like interactions. To high¬ 
light this, in the inset of Fig. 4 we plot the SE products 
DsTj and DgT*, where r* is the collective relaxation time 
calculated at the peak of S{Q), which is found to give 
the dominant contribution to viscosity [29]. Results are 
shown there from available experimental and numerical 
data points, as a function of the distance to the glass 
transition e; the latter being defined as e = 1 — 4)/4>g 
for experiments and as 1 — <?^'th/0th g simulations. 
Here, 0 th g = 0.235 ±0.010 and 0.215 ±0.010 are the ex¬ 


trapolated MCT glass transitions using a power-law fit 
for A^agg = 120 and 500, respectively [52]. While data 
suffer from statistical noise, it is clear from both sets of 
data that no systematic deviations from SE up to the 
largest studied volume fractions are observed. These re¬ 
sults cover a range which even exceeds the MCT tran¬ 
sition (e < 0) for the BD data, thus SE validity is pre¬ 
served even at very large degree of supercooling where 
deviations should normally be large. It is interesting to 
note that this behavior is consistent with previous obser¬ 
vations that star-polymer glasses, differently from other 
colloidal glasses, undergo a cessation of the aging process 
after a certain time 10^ s) [53] and in fact achieve equi¬ 
librium, which could provide a physical explanation of 
the suppression of dynamic heterogeneities in these sys¬ 
tems. Our findings for an ultrasoft coarse-grained model, 
which quantitatively reproduces the structural and dy¬ 
namical behavior of frozen block copolymer micelles as a 
realization of star-like soft colloids, are in agreement with 
those of simulations of the Gaussian core model, another 
effective potential often used to describe the behavior of 
soft particles [32] , but are in stark contrast with both ex¬ 
perimental and numerical results on hard-sphere colloids 
[27, 29]. Thus, it appears that softness is at least a neces¬ 
sary requirement in colloidal systems for the suppression 
of dynamical heterogeneities and the persistence of the 
validity of SE relation close to the glass transition. 

To summarize, we have examined the dynamics of soft 
colloids in the vicinity of the glass transition by a combi¬ 
nation of experimental techniques for different values of 
particle softness. Through BD simulations of a coarse¬ 
grained, ultrasoft model [40], we have been able to accu¬ 
rately describe the measured dynamic structure factors 
S{Q,t), both as a function of the micellar volume frac¬ 
tion and as a function of the wavevector, for different 
values of micellar aggregation number. We also observed 
that the increase of the macroscopic viscosity and of the 
inverse mesoscopic, self-diffusion coefficient is quantita¬ 
tively captured by the theoretical data. The coherent 
increase of both quantities indicates no violation of the 
SE relation up to the glass volume fraction. Our findings 
impressively confirm the microscopic origin of colloidal 
“softness” and its effects on the validity of the Stokes- 
Einstein relation for degrees of metastability for which it 
normally breaks down in the case of hard colloidal and 
molecular systems. In this way, they open up new realms 
for understanding and tailoring complex fluids not only 
with respect to structure and phase behavior [11] but also 
for colloidal dynamics [54]. 
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